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Abstract 

The operational characteristics of the Exhaust- 
Modulated Plasma Rocket are described. Four basic 
human and robotic mission scenarios to Mars are 
analyzed using numerical optimization techniques at 
variable specific impulse and constant power. The 
device is well-suited for “split-sprint” missions 
allowing fast, one-way low-payload human transits 
of 90 to 104 days, as well as slower, 180-day, high- 
payload robotic precursor flights. Abort capabilities, 
essential for human missions, are also explored. 

1. Introduction 

Notwithstanding the present emphasis on the 
assembly of a permanent space station in low Earth 
orbit, human missions to the planet Mars and beyond 
continue to captivate the imagination of scientists 
and engineers. Such missions embody the ultimate 
expression of human adventure and exploration. To 
be sure, steady, albeit limited, efforts to map out 
human exploration scenarios have continued. '• 2 - 3 
These scenarios have stimulated many dormant 
technologies and helped integrate them into 
increasingly realistic vehicle propulsion concepts. 

Two areas of special interest are space nuclear 
power and advanced (nonchemical) propulsion. 
Foreseeable developments in these areas will, 
nevertheless, still result in systems which are 
inherently power-limited. For example, in the 
development of advanced propulsion, the attainment 
of high exhaust-velocity or specific impulse (Isp) 
comes at the expense of vehicle thrust. Gains in 
rocket performance, measured as payload mass 
fraction, imply very long trip times. Conversely, 
short trip times result in very low-payload capability. 

This unfortunate situation can be greatly alleviated 
by the continuous modulation of the rocket exhaust, 
permitting considerably shorter transit times with 
reasonable payload. This technology has now come 
of age with the development of plasma heating and 
magnetic confinement schemes for fusion research. 
These combinations have been successfully adapted 
to rocket propulsion and have resulted in a high- 
power density rocket which is also capable of 
continuous exhaust modulation at constant power. 4 
The performance and capabilities of such a tunable 
plasma device are the subject of this paper. 

2. Operational Modes 

The concept of exhaust modulation has been known 
theoretically since the early 1950s. 5 6 ’ 7 Using this 
technique, maximum payload capability for a given 
trip time can be obtained through an optimum 


schedule that takes into account the local strength of 
the gravitational well in which the vehicle moves. 
Moreover, such systems are inherently flexible and 
provide the rocket with not only propulsion during 
the cruise phase, but also maneuvering capability 
upon reaching its destination. These systems also 
allow the implementation of a two-stage “split 
sprint” mission (see section 4) with the same rocket. 
One-way, 180-day nonhuman cargo missions with 
66% payload can be achieved. At the same time, 90- 
to 100-day fast human missions with from 2 to 14% 
payload capacity are also possible. 

Until the late 1960s, the technology to construct 
these systems had remained elusive. Nuclear-electric 
rockets suffering from limitations in power density 
due to electrode erosion focused on the achievement 
of high, fixed Isp. Chemical and nuclear-thermal 
rockets with fixed area nozzles, on the other hand, 
are incapable of thrust/lsp variation at constant 
power. In contrast, the present radio frequency (RF)- 
heated, magnetically vectored plasma design does 
not suffer from these limitations. In this type of 
rocket, exhaust parameters are mainly determined by 
input power and the confining properties of the 
magnetic field. 

The utilization of magnetic nozzles alleviates the 
materials constraints which impose limits on exhaust 
temperature. Such rockets exhibit a wide operational 
range in two basic modes: a high thrust/low Isp 
initial profile, followed by a low thrust/high Isp one 
as the vehicle loses mass and escapes the 
gravitational well. 

3. System Description 

As described elsewhere, 8 ’ 9 ’ 10 the Variable Isp Plasma 
Rocket is an electrodeless, electrothermal, RF- 
heated plasma propulsion device capable of 
continuous thrust/lsp modulation at constant power. 
The system, shown in Fig. 1, consists of three major 
magnetic cells, denoted as “forward,” “central,” and 
“aft,” all of which are linked together to form an 
asymmetric tandem mirror machine. 11 The forward 
end-cell comprises the main injection and ionization 
subsystem; the central-cell serves as a power 
amplifier. The aft end-cell provides the flow 
vectoring and exhaust and ensures its efficient 
detachment from the magnetic field. The general 
system parameters are listed in Table 1. 

Operationally, neutral gaseous fuel (typically 
hydrogen) is injected at the forward end-cell and 
ionized there by electron cyclotron resonance 
heating. This initial cold plasma is subsequently 
heated and controlled to the desired temperature and 
density by ion cyclotron resonance heating. The 
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Figure 1. Schematic representation of the variable Isp (or exhaust-modulated) plasma rocket 


latter process occurs in the much larger central-cell 
of the device. After heating, the plasma enters a 
two-stage hybrid nozzle at the aft end-cell where it is 
exhausted to provide modulated thrust. 


Table 1. General System Parameters 


Variable frequency 

500 kHz to 10 Mhz 

Four double half loop 
antennae in central-cell 

2 MW - each 

Two double half loop 
antennae in end-cells 

1 MW- each 

Transmitter operating 
voltage 

10 - 30 KV 

Central-cell magnetic field 

0.24 T 

End-cell magnetic field 

0.6 T 

Minor field 

1.2 T 

Mirror ratio 

2.0 

Superconductor 

NbTi or 

i high temp. S.C. 


The basic superconducting tandem mirror geometry 
has been modified to exploit the inherent axial 


asymmetry of these devices in order to provide a 
preferred flow direction. This approach is used to 
guide, control, and radially confine the plasma from 
injection to exhaust. The tandem mirror also ensures 
magnetohydrodynamic stable operation over a wide 
range of plasma temperatures and densities. 
Additionally, the mirror enables exhaust Isp and 
flow-rate control through manipulation of the end- 
cells as both input and output controlling gates. 

Experimental and theoretical investigations of this 
propulsion concept have been carried out over 
several years. 12 - 13 - 14 Based on these results, a 10MW 
nuclear electric system operating at 68% efficiency 
has been considered. In terms of specific power or 
“Alpha,” a value of 6 Kg/KW has been assessed. 15 
This value, although somewhat optimistic, reflects 
brisk progress in nuclear electric technology, as well 
as the high-voltage (low current) characteristic of 
RF-heated propulsion. The expected rocket 
performance parameters are shown in Fig. 2. 

4. Mission Concept 

The mission approach considered here is the familiar 
“split-sprint” scenario envisioned in numerous 
previous studies. A one-way, slow, high-payload 
capacity, automated cargo ship leaves first and 
places a habitat, fuel, and supplies on Mars. Some of 
the infrastructure will be located at strategic points 
on orbit and on the planet’s surface (an operating 
nuclear power plant on the surface, fuel and 
redundant landers on Mars orbit). Additionally, a 
fuel supply will be located on orbit at the edge of the 
Earth’s gravitational sphere of influence. 







































Figure 2. Theoretical rocket performance envelope for 
hydrogen propellant at 10 MWe and 60% efficiency. 


Following completion of this robotic mission, a 
smaller, low-payload, fast ship will carry the crew to 
Mars. The ship will pick up its final fuel load for 
Mars at a staging rendezvous point at the edge of the 
Earth’s sphere of influence. Failure to do so will 
abort the mission back to Earth. The fast cruise to 
Mars will be under modulated exhaust at maximum 
power on a trajectory designed for Mars-orbit 
insertion. Rendezvous with supplies and lander at 
Mars will be required to effect a landing. 

An important feature of a human mission should be 
its abort capability. Failure of a major system while 
en route could require return to Earth or, in some 
cases, use of the destination planet as a “safe haven.” 
Such abort capability was inherent (and used) in 
project Apollo and exists today in all Space Shuttle 
missions. Accordingly, interplanetary human 
missions should also possess such capability. 

Finally, as shown in section 6, the relative motion of 
Earth and Mars poses an interesting challenge to the 
mission designer. Very fast human opposition-class 
(opposition is defined here as the opposite alignment 
of the Sun and Mars with respect to an observer 
situated on Earth) missions are possible with exhaust 
modulation; however, a symmetric return leg cannot 
be accomplished in the same Martian year for transit 
times of less than 104 days. Symmetric outbound 
and inbound missions of 90 days (and possibly less) 
one-way are possible if the crew “winters down” on 
Mars for a period of 2 earth years. 

Such scenarios, while far too ambitious for the first 
mission, are not implausible for future expeditions 
building up to a permanent outpost on the planet. 

The use of Martian raw materials and installation of 


abundant surface nuclear power will more than 
enable these long stays. 

5. Mission Analysis 


The power-limited equations for the interplanetary 
vehicle will be described here briefly. More detailed 
descriptions are provided by Irving and Melbourne. 16 

The thrust, ‘T,” the power in the thrust beam, "p,” 
the specific mass of the power plant, “a,” and the 
thrust acceleration, “a r ” for power-limited systems 
are written. 


T = m p c 

1 • 2 
p = ym p -c 



a r 


_ ~m p c 
m 


(l) 



e 


where m p is the propellant flow rate and c is the 
exhaust velocity. This quantity is related to the 
specific impulse “Isp” through the familiar relation: 
c = g.Isp, where g is the acceleration of gravity at 
the Earth’s surface, 9.8 m/sec 2 . The remaining 
quantities are the total spacecraft mass, “m,” the 
power plant mass, “m r ” and its power rating, ‘7V’ 


A rocket equation for power-limited propulsion 
systems can be written from the above equations by 
solving for “(a r ) 2 ” and then integrating over time as 
follows: 


7 

m(t) 



( 2 ) 


where “m(rj” is the spacecraft mass at time, “m" 
is the initial spacecraft mass, and is the power 
plant efficiency. 


The vehicle parameters can be separated from the 
trajectory parameters and the rocket equation can be 
rewritten as: 


7 

m(t) 



m 0 m w 


where 


P(r ) 2 


a 

2~e' 



■dt 


(3) 
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At the final time, “t f , ” “p” can be written as, 


W r /)] 2 = f y2 ’ where: 

J = 2 

is the performance index. 

The trajectory is optimized independently of the 
vehicle by minimizing J 2 , which effectively 
minimizes the amount of propellant required for the 
mission. This optimized trajectory is then used with 
the vehicle parameters, which are determined by the 
mission planner, to calculate the power plant and 
payload mass fractions. 

The spacecraft mass at any time, “t, ” is equal to the 
sum of the payload mass, “m t ,” the power plant 
mass, “m r ” and the propellant mass, “m (t). ” 

m(t) = m L + m w + m p (t) (4) 

The initial propellant mass, “m (0), ” is assumed to 
be exactly the amount of propellant necessary to 
complete the mission, so at the end of the mission, 
m (t.) = 0. 

From reference 16, the maximum payload mass 
fraction and the propellant mass fraction can be 
shown to be. 


m, 


m r 


= (1 "tf 


(5) 


Isp(t) = 


1 

g-%(t) 



+ J • dt 

^0 


( 8 ) 


To compute the optimal Isp as a function of time, the 
optimal control functions "a Tx> ” “a Ty ," u a Te ” and the 
resulting value for “7” need to be found. The 
functions “fl K ” “a T} , ” and are found by 
solving an optimization problem defined as follows. 

The performance function is, 

J = — • tij-) • dt (9) 

*0 


The equations of motion are: 
r = v 

( 10 ) 



where |i is the Sun’s gravitational parameter, 1.32715 10 20 
m 3 /sec 2 . 

The initial conditions are: 


?0 = f e( t oX 

v 0 = v e (t 0 ) 


( 11 ) 


where r e and v e are the position and velocity vectors 
of the Earth at the initial time, t 0 . The final 


m p (Q) 


m 


o 


=P 


( 6 ) 


=r m (t f ) 

v f = v m (t f ) 


( 12 ) 


where 


P= 



conditions are: 

where ? m and are the position and velocity 
vectors of Mars at time, t f . 

From the calculus of variations, as shown in Bryson 
and Ho, 17 the Hamiltonian is defined as 


Another parameter of interest is the Isp time history 
of the optimized transfer. Isp(r) is a function of the 
thrust acceleration time history, the integral of the 
thrust acceleration squared, a, e, and 7. 


H = i>,-fi + L (13) 

i=l 
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where L = ^a T 2 from equation (9), “/’’are the right- 

hand side of equation (10), and “A ” are the Lagrange 
multipliers corresponding to equation (10). 

From the first necessary conditions for an extremal. 



3x(X-r) ^ 







(14) 


where, 

x.=MA) 


and 



(15) 


After application of the calculus of variations, the 
optimization problem requires minimizing “J 2 " as a 
function of the initial values of the Lagrange 
multipliers, 

M0XM°)’M°) 

MoxVoxMO) (1< 


the final time, and satisfying the terminal 
constraints in equation (12). 

This constrained minimization problem is solved using 
a numerical optimization method in the program 
Hybrid Optimization Technique (HOT). HOT provides 
the inteiplanetary simulation by integrating equations 
(10) and (14), numerically using equation (15) for the 
control values, and using equations (11) and (16) for 
the initial state and multiplier conditions, respectively. 
The initial X's are iterated upon using the Davidon- 
Fletcher-Powell Penalty Function Method, 18 which also 
resides in HOT. 


6. Performance 

A representative set of mission scenarios, as 
described in section 4, was studied using the 
numerical techniques outlined above. No attempt has 
been made in this initial study to account for the 
gravitational effects of the departure and destination 
planets. Rather, the ships move in heliocentric space 
only under the gravity of the Sun. Boundary 
conditions require zero relative velocities between 
the ship and the planet at both the departure and 
arrival points. Furthermore, these points are chosen 
to coincide with the position of the appropriate 
planet at the departure and arrival times. Optimized 
trajectories accounting for the planets’ gravitational 
fields and low orbit departure and arrival, as well as 
assorted techniques for aerobraking and gravity 
assist will be taken up in a future study. 

Given the launch date and the desired length of the trip, 
the initial state vector of Earth and the final state vector 
of Mars (or vice versa) are calculated. From these 
boundary conditions, the optimization program 
computes the value of which, roughly described, is 
a measure of the fiiel necessary to accomplish the 
given mission. From the above formulation, the 
assumed values of “a” and “e” are used to calculate 
the value of “(3,” the propellant mass fraction of the 
rocket. From “(3,” the payload and power plant mass 
fractions are also calculated. The trajectory is 
computed considering only the heliocentric sphere of 
influence, and seeks to match perfectly the initial and 
final state vectors of the planets. 

A 90-day trajectory profile from Earth to Mars (Fig. 
3) will be used to illustrate the performance and 
properties of an exhaust-modulated rocket. The 
trajectories are integrated in convenient half-day 
(43,200 seconds) segments. Vehicle moves in 
heliocentric space matching the initial and final 
positions and velocities of Earth and Mars, respec¬ 
tively. Vehicle starts with zero Earth-relative 
velocity and arrives at Mars with zero relative 
velocity as well. Techniques to optimize Earth 
escape and Mars capture are considered separately. 
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Figure 3. Typical one-way 90>day Mars transit with 
exhaust modulation. 

The thrust acceleration vectors associated with the 
90-day trajectory are shown in Fig. 4, while the Isp 
as a function of time is plotted in Fig. 5. Inspection 
of these graphs shows the greatest velocity changes 
near the beginning and the end of the mission. The 
midcourse phase of the trajectory is a high Isp, low- 
thrust “quasicoast” period. The vehicle mass 
variation as a function of time is shown in Fig. 6. It 
indicates the greatest fuel expenditure at the 
beginning and at the end of the trajectory. 



the radial distance, the more the rocket must struggle 
against the gravity well of the Sun. Therefore, the 
best possible one-way trajectories arrive at Mars, or 
are launched from Mars, when the planet itself is 
near perihelion. This is shown in Fig. 7. 





Figures. Isp profile over tbe entire transit Isp and thrust 
are exchanged continuously at constant (maximum) power. 
Hie profile is nearly symmetric showing only a minor 
difference at Mars, due to the combined effects of a lower 
Sun’s gravity and reduced total mass at arrival. 



spent during the initial and final high-thrust phases of 
the mission. The use of gravity assist and aerobraking 
techniques may further reduce the fuel expenditure at 
these two ends. 


Figure 4. Trajectory plot for a 90-day transit showing 
the thrust acceleration vectors at 10 intervals during 
the mission. Thrust acceleration is nearly radial, 
decreasing in magnitude toward the midpoint of the 
mission, reversing direction, and increasing again at 
Mars arrival. 

The eccentricity of Mars’ orbit (.0934) causes a 
considerable variability in the radial distance an 
Earth-Mars trajectory must cover. The Earth’s orbit, 
being much less eccentric (.0167), is less of a factor 
when determining optimal trajectories. The greater 


One can define a phase angle i) for the transit as the 
angle subtended between two radii originating at the 
Sun and intercepting the departure and arrival points. 
Trajectories which traverse small phase angles are 
generally fuel-inefficient, as confirmed by experi¬ 
mentation with these approaches. A phase angle of 
roughly 60° yields the optimum results for a 90-day 
mission. Phase angles more or less than the optimum 
values, as determined through iterative runs, are 
fuel-inefficient. Thrust is wasted to counteract 
components of the initial velocity which would 
improve, rather than hinder, the rocket’s travel time 
to Mars. 
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Figure 7. The eccentricity of Mars’ orbit provides an 
optimum window of opportunity for the transit when 
the planet’s orbit is closest to that of Earth. 

The Martian year of 686.9 days forces a period of 
opposition of the planets approximately once every 2 
years. In the years 2016 and 2018, the oppositions 
are relatively close to the perihelion of Mars. These 
2 years offer three launching options: two human 
and one robotic. For example, the “tugboat” profile 
shown in Fig. 8 gives little regard to the travel time 
of the rocket and seeks to maximize the payload 
mass fraction placed in Mars orbit. As described in 
section 4, such a profile is used to deliver the habitat, 
return fuel, and supplies to Mars. Once on the planet, 
in situ resource utilization machinery could be set in 
motion. The nuclear power plant on board might 
also be used to provide supplemental power to the 
habitat and, hence, become part of the useful 
payload upon arrival. 

Once the operation of the equipment is started and 
verified from the safety of Earth, the human mission 
could commence during the next 2-year cycle. The 
profile shown in Fig. 8 is a 180-day mission with a 
66.66% payload mass fraction. However, longer 
missions with better mass fractions are possible. 

The first human mission considered here is a short 
one shown in Fig. 9. This profile is executed within 
the same Martian year. The mission encompasses a 
101-day outbound trip, a 30-day stay, and a 104-day 
return This “speedboat” option is one of the shortest 
Mars missions thus far projected, albeit with a 2% 
payload mass fraction in each direction. 

Such short human transits ameliorate the 
physiological hazards associated with long human 
exposure to the space environment. Moreover, the 
continuous acceleration provides an artificial gravity 


of sorts which may have added positive effects. 
Reduced mission duration (235 days) would put 
much less of a requirement on the development of 
long-duration facilities and their associated 
maintenance. Radiation exposure would be reduced 
as well. 


Trajectory Plot 



Figure 8. Typical trajectory plot for a high payload, 
longer duration transit Values of 18% payload-mass 
fraction and 90 days are also possible with the same 
engine, illustrating the dual nature of this propulsion 
system. 


Trajectory PJot 



inbound and outbound legs. While shorter transits are 
allowed by the propulsion system, the planet’s 
alignment becomes unfavorable for transits less than 
100 days. Shorter transits are possible if the mission 
includes an extended stay on the surface of the planet. 

Reducing the time of stay on Mars in this profile 
does not yield substantial improvements in the 
transit times or payload mass fractions. Shortening 
the stay to 20 days might yield a decrease of 6 days 
outbound and 5 days inbound and keep the payload 
mass fraction constant. 
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7. Abort Scenarios 

An important attribute of the exhaust-modulated 
plasma rocket is its powered abort capability. 

Should the habitat fail, the crew become 
incapacitated, or one or more of the modules of the 
rocket cluster fail, the system has a limited 
capability to return directly to Earth or abort to 
Mars, albeit with an increased travel time. 

The issues associated with powered aborts back to 
Earth are the nominal mission profile, the time of 
abort, and the abort duration. Given the nominal 
mission profile and the time at which the decision to 
abort is implemented, the fuel remaining and the 
initial state vector of the spacecraft are determined. 
The abort duration specifies the location of Earth at 
the arrival point and, consequently, the final state 
vector. The same optimization program is then used 
to determine the necessary mass fractions to 
accomplish such an abort. Comparing the mass 
fractions of the vehicle at the time of abort and the 
mass fractions necessary to accomplish the abort, it 
is then possible to determine whether the desired 
abort is feasible. 



Figure 10. Typical powered abort returning the 
spacecraft to Earth 75 days after abort condition is 
declared. Such abort windows exist within 15 to 20 days 
after Earth escape. The scenario depicted here assumes 
a full-up propulsion system after the abort is declared. 

The abort scenarios considered here are shown in 
Figs. 10 and 11. They occur 15 days into the 90-day 
mission, and last 75 and 180 days, respectively. Both 
are possible; however, the 180-day return can make 
it back with fuel to spare, as in a propellant system 
failure. Conversely, if the fuel is available, this 
demonstrates the possibility of aborts even after 15 


days en route. In these cases, the velocity away from 
Earth orbit and the fuel remaining are the driving 
factors. In general, the efficiency of fuel usage 
increases with the abort duration. While the 
calculations suggest the possibility of powered 
aborts due to a partial engine cluster failure, these 
have not yet been simulated. 

The three scenarios depicted above embody the 
fundamental operational requirements for a human 
mission to the planet Mars. 



longer return patch would result from a propellant 
system failure reducing the fuel available and hence 
increasing the return time. Abort windows exist within 
15 to 20 days from Earth escape. 


JD-'24$B287,00 day* outbound. JO.24S90M 5,90 dtya ratum 



symmetric inbound and outbound legs. Payload-mass 
fractions are considerably higher with this option, due 
to the optimum alignment of the planets. 
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8. Extended-Stay Missions 


10. References 


Finally, the extended-stay Mars mission is 
considered. The complete profile is shown in Fig. 

12. The 90-day outbound trip is launched during the 
optimum point in an initial Martian year. This choice 
increases the payload mass to 18.12%. After a 
704.5-day stay, the craft returns in 90 days with a 
14.49% payload-mass fraction. The short transit 
times minimize radiation exposure in space, 
assuming improved shielding of the Martian habitat. 
The improved payload-mass fractions are possible 
because of the 2-year wait for the next opposition of 
Earth and Mars near Mars’ perihelion. 

Such a long-stay mission would require refinements in 
closed-environment recycling techniques, supple¬ 
mented by extensive in situ resource utilization. How¬ 
ever, the benefits gained are payload mass fractions 7 
to 9 times better than those offered in the 30-day stay 
mission profile. Perhaps this would be the choice 
profile for a second, longer duration mission to Mars, 
once the lessons from the first have been absorbed. 

The trajectories examined have not accounted for 
potential improvements in performance, such as 
aerobraking, gravity assist, and capture into a highly 
elliptical orbit. Currently, use of the “crew taxi” 
concept is assumed. The crew would be flown out to 
the main vehicle once it had approached escape 
velocity. However, given its high-thrust capability, 
further consideration of the performance of the 
rocket in Earth’s and Mars’ spheres of influence is 
planned. 

9. Conclusion 

For power-limited rockets, the technique of exhaust 
modulation at constant power provides the optimum 
approach to fast interplanetary missions. A “split- 
sprint” mission scenario can be accomplished with a 
one-way, 180-day robotic space-tug precursor 
mission, followed by a 101-day, low-payload, 
human fast boat. After a 30-day stay, the human ship 
can return to Earth in 104 days. For a long-stay 
mission, the transit times can be considerably 
shorter and the payload-mass fraction increased. 

This modality may be utilized when a permanent 
habitat is in place. 

Such exhaust-modulated plasma rockets have now 
reached the experimental stage and have become 
increasingly attractive as a result of recent 
developments in ion cyclotron plasma heating 
techniques and superconducting magnetic confinement 
and vectoring systems. Furthermore, these engines 
provide high-thrust orbital maneuvering capability and 
abort margins which are essential for human missions 
to the planet Mars and beyond. 
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